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1 Introduction 



The aim of this paper is the computation, within algebraic Bethe ansatz method, 
of the spontaneous magnetization of the XXZ spin-i chain (see Q>@])> ^ taking the 
thermodynamic limit of the form factor formulas of local spin operators for the finite 
chain obtained in ||. 



We consider the XXZ Heisenberg chain of finite length M 
which is given by 



], the Hamiltonian of 



M 

H A = { ff mCi + ^mCl + A ( CJ m cr m+l ~ 



(1.1) 



m=l 



with periodic boundary conditions 



a i — a M+li 



x,y,z. 



This Hamiltonian acts in the quantum space TL of the chain, which is the tensor product 



of M local quantum spin-i spaces 7i m isomorphic to C 2 and a° 



x, y, z, are the 



standard Pauli matrices acting in TL m . In the anisotropic case (A ^ ±1) this model was 
solved by means of the Bethe ansatz in ||[?J (see also |§). 

We use the algebraic Bethe ansatz solution for this model II, taking the XXZ R- 
matrix in the following normalization (as a matrix acting in C 2 <8> C 2 ): 



( 1 \ 

6(A,/x) c(A,/u) 

c(A,/i) b{\,n) 

\ 1 / 



(1.2) 



The functions b(X, fi) and c(A, fj,) are defined as 

sin(A — fi) 



b(X,ti) 



i sinh (, 



sin(A — ii — iQ ' ' sin(A — (J, — iQ ' 

where the parameter £ is related to the anisotropy parameter A of the Hamiltonian as 

A = -(q + q ), withq = e^. 

The i?-matrix is a linear operator in the tensor product of two two-dimensional linear 
spaces V\ ® V2, where each Vi is isomorphic to C 2 , and depends generically on two 
spectral parameters Ai and A2 associated to these two vector spaces. It is denoted by 
-Ri2(Ai, A2). 
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The monodromy matrix is constructed as an ordered product, in an auxiliary space 
Vo isomorphic to C 2 , of such i?-matrices Ro m acting in Vo TC m : 

T(\)=R OM (\ + i^...R 02 (\ + i^)R 01 (\ + i^) = ( £M 5(A) 

In the last formula, the monodromy matrix is represented as a 2 x 2 matrix in the 
auxiliary space Vo, whose entries A(\), B(X), C(A), and D(\) are operators in the 
quantum space Ti of the chain. 

Since we will further use some results of ||, let us note, to compare notations, that 
the spectral parameters are equal here to the corresponding spectral parameters of || 
multiplied by i, and that £ corresponds to —rj of ||. Moreover, we consider here the 
homogeneous case where all the inhomogeneity parameters £ m of §§ are equal to — (/2. 



The eigenstates of the Hamiltonian (|1.1|) can be constructed by the action of the 
operators B(\) on the ferromagnetic state | 0) (which is the state with all the spins up). 
More precisely, the state B(Xi) . . . B(\n)\ 0) is an eigenstate of the Hamiltonian ( |1.1[ ) 
if the set of spectral parameters {^j}i<j<N is a solution of the Bethe equations 

f M^ + f Yfr MXj-h-iQ x< (13) 

y sin{Xj -£)) l} lS m(X j -\ k + ^0 ~ J - 
The norm of Bethe eigenstates is given by the Gaudin formula |8],[l0| (see also Q): 

(0| n^(A,) n^(A*)|0> = (-sinhC)^ II — (x \ , det.V({A a }), 

3=1 fc=l a^fS Sm[Aa A(3) (1.4) 

where M is an N x iV matrix the elements of which are given by 

.» , f f sm(A„ - f ) \ M A gjgtj. lA+ '0 1 

^=-W n {U^7fyj n sin(Aa _ At _ iC) }- a-*) 

The ground state of the XXZ model in the region A = ^g + g -1 ) > 1, q > 1, is 
degenerated in the thermodynamic limit (M — > oo), namely there are two states with 
the same energy which we will call the ground state | ^i ) and the quasi-ground state 
| ^2 ) (° n the finite lattice, these states possess different energies). These two states can 
be obtained in the algebraic Bethe ansatz framework (see 11). 

The spontaneous magnetization is defined as the modulus of the normalized matrix 
element of the local spin operator between these two states: 



so 



(*lKI*5 



i , _ i _ , I 

2 



(#1 1 #1)2 (*2 I ^2) 



(1.6) 
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In p|,|Tl||, Baxter proved that the spontaneous magnetization of the infinite chain could 
be expressed as an infinite product: 



so 



JJlT^J , (M = oo). (1.7) 
This formula in the thermodynamic limit was also reproduced by means of the q-vertex 



operator approach (see |12| ). 

In this paper, we study the spontaneous magnetization so in the framework of the 
algebraic Bethe ansatz. From the results of paper || , the corresponding form factor on 
the finite chain is represented as a ratio of two determinants of M x M matrices. In 
the thermodynamic limit, it becomes the ratio of two Fredholm determinants of linear 
integral operators, which can be computed explicitly. This leads directly to the Baxter 
formula (jOj). Thus, we obtain here an independent proof of the Baxter formula from 
algebraic Bethe ansatz. We prove also that there are no perturbative (j^hz) finite size 
corrections so that for the finite chain the Baxter formula has only exponentially small 
corrections. 

This paper is organized as follows. 

In section |2[ the expression for the form factor of operator obtained in Q is 
recalled. It is reformulated in a form suitable for the thermodynamic limit and partic- 
ularized to the ground and quasi-ground states. 

In section ^, the procedure of taking the thermodynamic limit is discussed. 

In section |4], the spontaneous magnetization is represented in terms of Fredholm 
determinants, which gives a proof of the Baxter formula in the thermodynamic limit. 
The finite size corrections are also evaluated. 

2 The form factor 

Let us consider the function 

N N 

<oi n c( Mi x n B(\ k )\o) 

°*W = ^ > (2- 1 ) 

(0|riC(A,) UB(X k )\0) 

3=1 k=l 

where {Afc}i<fc<Ar and {^j}i<j<N are solutions of the Bethe equations 

11 sinQ, — \, 4- it\ ~ V ' 



v sin(Aj - f )J lJ[ sin(Aj - X k + i() 
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These equations can also be written in a logarithmic form: 

N 



M Po (X,) + d ( X i ~ X k) = 2 ™i> 1 <3<N, (2.3) 



fc=i 



where rij are integers for N odd and half integers for iV even. The bare momentum 
Po(X) and the scattering phase 6(X) are defined as 



sin(A + f 

6{X) = iln 



Po (X) = i\n 

sm(A 1 1 



2 ) 

sin(z£ — A) 
sin(i£ + A) 



In the recent paper ||, an expression for any arbitrary matrix element (form factor) 
of the cj^j operator between two Bethe states has been obtained in the framework of 
algebraic Bethe ansatz. Using this result (see proposition 5.2 of ||) as well as the 
Gaudin-Korepin |8, 10] formula for the norm of Bethe states, one can represent the 
function s z (m) as follows: 



|i sin(/z fc - Hj) detM '^^ 



./ -nX— /> \ , NN ] rr sin(At. — A,) det(<S — 2Q) 
»-jno- cxpi -i{m-l) 2^ (po(Afc) -Pofjik)) \ [[ 

h=\ 

where the N x N matrix S is given by 

$ ± = sinhC t^i- sin(A b - \x k ± i() 

ab sm(X b - Ha) sin(A fe - \i a ± i() ^1 s in(A fe - X k ± <) ' 

and Q is the following N x N rank one matrix: 

sinhC J^sin(/i fc -i|) 
Qab = ~7T~ ~lc^~77. JcTll 

A/" is the Gaudin matrix: 



sin(^ a + i|) sin(^ a - i|) ^ sin(A fc « 2 



iV 



AA afc = 5 a6 {M^(A a ) - 2vr^^(A a - X k )} + 27tK{X a - X b ), (2.5) 



k=l 



4 



where the function K(\) is proportional to the derivative of the scattering phase 

K( X) = ->(A) = - ...f?/* ... . (2.6) 
Z7T 2tt sm(A + iQ) sin (A — iQ 



Note that the fact that Q is a rank one matrix together with the orthogonality of two dif- 

sin(Mfc-»C/2) 
sin(A fc -iC/2) 



ferent Bethe states (i.e. det 5 = 0) allowed us to insert the global factor Y\ k=1 ^™^ fc ^/ 2 ) 



in the matrix elements of Q. 

Thus, we will derive the Baxter formula from the expression ( |2.4| ). In order to 
take the thermodynamic limit in (|2.4| ), we first rewrite it in a more convenient form by 
introducing the following matrix 

TV 

ef(Ao _ w) UMvb-Xk) 

Xab = "T— 77 r i-r . 7 c> (2-7) 

sm(A a -fi b ) [[ sm((j, b - n k ) 

k^b 



the determinant of which is 



det X = exp < i ^^(Aj — fij) ? 



N EI sin(A fc - Xj) 

j>k 



J] sin(^ fc - fJ-j) 
j>k 



Hence ( |2.4| ) can be rewritten as 

N N f S 

s z (m) = exp < — i(m 



!)X]Oo(A/c) -Po0*fc))| exp ~ Aj)| —^^7 ^ 

fc=l ^ ^ 7=1 ^ 



The products of matrices M = A'iS and V = XQ can be calculated using some identities 
for the rational functions. The procedure is similar to the one used in Appendix A. Here 
we merely give the results: 

M ab = -2itK{\ a - At) + 

II sin(A a -A fc ) , N N . x 

fc^a / -r-r sm(A a - Hk + iQ i-r sm(A a - fi k - iQ \ 



n S m(A a -^) \L\ shl ^-^ +i o lAMXa-\ k -iOj ' (2 . 

fc=l 

./exp{i(A a -i§)} exp{i(A a + 4)} /-vv , n 
V " = \ sin(A a -4)) " sm^ + 4)) ^ t'g^ " ^ , 



(2.9) 
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Note that these expressions are valid for any two different solutions {Aj}i<j<7v and 
{^fc}i<A,<Af of Bethe equations. To calculate the spontaneous magnetization ( |1.6| ), one 
has to particularize them to the ground and quasi-ground states | } and | ^2 )> that 
is to compute the quantity 

- (m) = (»/«!,) ■ (2 - 10) 

From now on, {/xi, . . . /j>n} and {Ai, . . . \n} will denote the sets of spectral param- 
eters solutions of the Bethe equations (^) corresponding respectively to the ground 
state I ^1 ) and the quasi-ground state | ^2 )• For both states, N is equal to M/2. The 
ground state | ^1 ) is parametrized by the following set of nf 

N + 1 

nj = — +j, j = l,2,...,N, (2.11) 

whereas the quasi-ground state | ^2 ) is parametrized by the shifted set: 

N+ 1 

h j = j = l,2,...,N, (2.12) 

(see, e.g., We will moreover use the following notations: 

5Xj = Xj+i — Aj, 5/j,j = Hj+i — (Aj , <5Aj = fij — Xj . 

The energies of these two states are equal in the thermodynamic limit, the difference of 
their total momenta being equal to tt. 

Thus, one has expressed the quantity fl2.10|) in the form: 

..(m) = (-ir-^ P \i±fX, y et ^- r 2V \ (2.13) 



detTV 



where A4 and Af are the N x N matrices the elements of which are given by (2.8) and 
(|2.5|) , whereas the elements of the matrix V are simply (since the difference of the total 
momenta of the ground and quasi-ground state is tt): 



e 2i\ a _ CQS fa £ 

sin(A a + «|) sin(A a - i- 2 , 

The expression ( p. 13 ) turns out to be very convenient to take the thermodynamic limit 
since the functional form of the non-diagonal terms of the matrices A4 and Af, as well as 
the matrix elements of V, do not depend on the number of sites M = 2N of the chain. 
Nevertheless, the diagonal terms of the matrices A4 and Af are more complicated and 
should be treated separately. 
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3 The thermodynamic limit 



3.1 General procedure 

In the previous section, we obtained an expression of the spontaneous magnetization 
as a quotient of two determinants. Most of the matrix elements involved are independent 
of N, which makes the thermodonamic limit (M = 2N — ► oo) quite obvious. The only 
terms we have to take care of (i.e. which depend explicitly of N) are on the one 
hand the global coefficient, and on the other hand the diagonal matrix elements of M 
and TV. Both kinds of terms consist of sums (or products) on the different spectral 
parameters in the ground or in the quasi- ground states. Thus, the first step in taking 
the thermodynamic limit is to know how to deal with such sums. This leads us to 
formulate the two following propositions. 

As previously, {^j}i<j<N and {^j}i<j<N denote respectively the set of spectral 
parameters corresponding to the quasi-ground state, solution of (|2.3|)-( |2?P2|) , and the 
set of spectral parameters corresponding to the ground state, solution of fl2.3|) -( |2.11| ). 

Proposition 3.1. Let f be a C°° ir-periodic function. Then the sum of all the values 
f(Xj) where the set of the spectral parameters {Xj}\<j<N parametrizes the quasi-ground 
state can be replaced by an integral in the thermodynamic limit according to the following 
rule: 

1 N f~ 

3=1 J ~% 



where p(X) is defined as the solution of the Lieb equation [lij 



,>(\) !- f ^ duK(X - v)p(y) = (3.2) 

tt/2 ^ 



The same equations (fO)-(3.2) stand concerning the sum on the values f(fJ-j) where 



{pj}i<j<N parametrizes the ground state. 

The proof of the proposition is given in the Appendix B. 

There exists a generalization of this result in the case that / is not itself a periodic 
function but that its derivative is periodic: 

Proposition 3.2. Let f be a C°° function such that f is ir-periodic. Then the sums 
of the values of f at the points Xj, 1 < j < N parametrizing the quasi-ground state, 
and respectively at the points fj,j, 1 < j < N parametrizing the ground state, are given 
in the thermodynamic limit as the corresponding integrals, the only rational finite-size 
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correction being of order tt : 



^E/( A i) = rd\f(\)p(\) + ^P- + 0(M-°°), (3.3) 
mJ2^) = j[dpf(p)p(p) + ^>+0(M~^), (3.4) 



3=1 

where c\ and c% are constants independent of M . 
3.2 The norm 



Using the proposition 3.1 and the Lieb equation (|3.2|) one can immediately calculate 



the diagonal terms of the Gaudin matrix J\f: 

Maa = 2ttM{p(A ) + jjK(P)} + 0{M~°°). (3.5) 
Thus the matrix elements of M have the following form in the thermodynamic limit: 
M ab = 2vrM{<W(Aa) + ~^ K ^a - Xb)} + 0(M-°°). (3.6) 

3.3 Diagonal matrix elements of M. 

In order to calculate the diagonal elements of the matrix A4, let us consider the 
following 7r-periodic functions: 

<P±(A a ) = I I — 7T — — r , (3.7) 

^ sm(A a - X k ± iQ 

N 

]1 sin(A a - p k ) 
4>{\ a ) = M y^ . — , (3.8) 

k^a 

and those defined similarly for the second set of parameters: 

A sin (/^ - A a ± «C) to n^ 

V>±(w0 = M — ( ^rpr, (3.9) 

sm{p k - p a ± iQ 

TV 

EI sin(// fc - A a ) 

^(/z*) = M^— (3.10) 

11 sm(p k - p a ) 

a^k 
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The diagonal elements of the matrix A4 (2.8) are actually: 

M aa = -2ttK(0) + iM<t>- l {\ a ) (<^+(A ) - <t>-{X a )) . (3.11) 
These functions are related by remarkable identities: 

1 sinh C , , , « i / , 

M sm(/i fc - A ffl ± «g sin(Mfe ~ -M 

the detailed proof of which, based on identities for rational functions, is given in Ap- 
pendix A. 

It is to mention that the factors 0(A) and VKaO were used earlier by N. Slavnov (we 
thank N. Slavnov for communicating us the results contained in his Ph. D. thesis). 

Summing up the identities ( |3.12 ) with plus and minus signs one obtains the following 
equation: 

2tt N 

— K (»k - A a )0(A a ) = tty+Vfc) - </C V*)), (3.13) 

a=l 



which can be rewritten in an integral form by means of proposition 3.1 



2 K(n - X)0(X)p(X)dX = ^+\p) - ^Z l {p)) + 0(M-°°). (3.14) 

2 

An analogous equation can be derived similarly for the functions ip(p) and <f>±(X): 

7T 

2 K(A - MMpiridn = ^-(4>+H X ) - </C X (A)) + 0(M-°°). (3.15) 

2 

The computation of the main order of the functions 0±(A), ip±(p) in the thermo- 
dynamic limit is rather simple but it is more complicated to prove that there are no 
rational perturbative finite-size corrections. We present this calculation in the Appendix 
C. The result is 

0±(A) = <M A ) = ±i + 0(M-°°). (3.16) 



Equation ( |3.14|) is used to calculate the function </>(A). By means of ( 3.16 ), it can 
be rewritten as 

K{jt - X)(f)(X)p(X)dX = -- + 0(M-°°), 



7T 
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which admits a unique solution 



0(\)pW = -- + o(m-™). 

7T 



Thus, one obtains the following expression for the diagonal terms ( [3.11 ) of the matrix 
A4 in the thermodynamic limit: 

M aa = 27rp(A )M - 2^fT(0) + 0(M-°°). (3.17) 

Hence the elements of A4 have the following form: 

Mob = 2nM{6 abP (X a ) - ^K{\ a - X b )} + 0(ikr°°). (3.18) 
4 Fredholm determinant representation 

In the thermodynamic limit, the determinants of the matrices M and Ai — 2V are 
replaced by Fredholm determinants: 

N 



detM= (2nM) N (flpiXj)) {det(I + K) + 0(M-°°)} , (4.1 
det(X-2V) = (2irM) N (]Jp{\ J )){det(i-k-2V) + 0{M~ 00 )}, 

V i= i / 



(4.2) 

where I is the identity operator and K is the integral operator acting on the interval 
[—§,§] with the kernel K(X — /x) (|2.6|). The kernel of the integral operator V is 

V(X,p) ^ V(X) = - : f ^ ~ ^° sh C (4.3) 
27r sm(A + «|J sm(A — &|J 

To calculate the coefficient in the representation ( 2.13[) one has to use the proposition 



3.2, and the first order of the Taylor developpement of 5Xj which has been computed 
in Appendix C: 

JV 



exp i i ^ S^j \ =i + 0(M' 

^ .7=1 ' 



3 = 

This leads to the following Fredholm determinant representation for the function 
s z {m) ( EIOD : 



s z (m) = i{-lT^ ^-J^ ^ + O(M-). (4.4) 
det(I + K) 
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The function 



s z (m) 



(* 2 K|*i) 



<*i|*i> 

can be computed similarly, which gives 

s z (m) = -s z (m). 

Therefore the spontaneous staggered magnetization admits the following represen- 
tation in terms of Fredholm determinants: 



•so 



det(J -K-2V) 
det(I + K) 



+ 0(M~ 



(4.5) 



What remains to do now is to compute explicitly these determinants in order to 
obtain the Baxter formula. This can be done by means of Fourier transformation. 

As the kernel of the integral operator K depends only on the difference of the two 
variables, and since we are in the regime q > 1, this operator can be diagonalized by 
Fourier transformation. Its eigenvalues are thus obtained as the Fourier coefficients k n 
of the function K(\): 

7T 

k n = [ 2 K(X)e 2iXn dX = e~ 2 ^ = q' 2 ^. (4.6) 

Therefore the determinant of the operator I+K is only the infinite (convergent) product 
of its eigenvalues, namely (q > 1): 

oo oo 

det(J + £)= [J (l + ^ 2|n| ) = 2[](l + ^ 2n ) 2 . (4.7) 

n=— oo n=l 

After Fourier transformation, the operator V becomes an infinite matrix which ad- 
mits only one nonzero column v n Q\ 

7T 

t>nO= I ' V(X)e 2iXn dX. (4.8) 



Since the operator I — K is diagonalized by the Fourier transformation, only the diagonal 
element vqq contributes to the determinant, namely 

oo 

det(7 — K — 2V) = -2v 00 [] (1 - q- 2n f . (4.9) 

n=l 
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The computation of the integral ( |4.8|) leads to vqq = — 1. Finally, from the expres- 
sion ( |4.5| ) for the spontaneous magnetization, one obtains the following result: 

(oo -. _2 n \ 2 
which coincides with the Baxter formula. 



Appendix A 

In this appendix we give the proof of identity ( 3.12j ) , which provides a typical example 
of a proof for identities involving rational functions. In particular, (2.8) and ( 2.14j ) can 
be proved similarly. 

We thus consider the following identity 



N 



sinh£ 



n sin(A a - fi r 

m=l 



N 

y 

^ sin(^ - A a + i() sm(fik ~ K) II sin(A a - A r 



N 

n 



sin(^ fc - n m + iQ 
j sin(/i fc - A m + <V A 1 



To prove this relation it is convenient to introduce the exponential variables: 

v k = e^ k , u a = e iXa , q = e ( . 
The left hand side of (A.l) can be represented using these variables as 



N 



m=l 



N 



1 



N 

n (ul - V. 



■ * V * 11 „ (^2 _ U 2 ){v 2 _ ^ U{U 2- U 



(A.2) 



For the right hand side one obtains 



r.h.s. 



N 

n 



j v m (yl - u 2 m q 2 



(A.3) 



Consider the sum in ( A.2 ) as a function of vj. 



N 



Wi(v 2 k ) = Y J 



q 2 -l 



N 

n («2 - « 

m=I 



^ K ~ U a)(^ - «a9 2 ) EI ( u a " «m) ' 
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which is to compare with the following function 



W 2 {vl) 



N 

n 



v 2 m q 2 



A 

m=l K 



U; 



2 n 2 ' 



oo one can easily see that Wj(v?) 



Both of them are rational functions. When v\ 

for J = 1,2. Also both functions have only simple poles in the points v\ = u^q 2 and 
their residues in these points are equal. Hence one can conclude that 

W x (yl) = W 2 (yH 



which leads immediately to the identity (A.l). 



Appendix B 



We present here the proof of the propositions 3.1 and 3.2. 

Let us begin with the proof of proposition 3.1. Note first that, in case of variables 
Qji 1 < j < A which are homogeneously distributed inside an interval of lenght it (i.e. 
qj+i — qj = it is easy to prove the analogous proposition (with p(X) = Jr): 

if.,. 1 



M ^ 



f{qj) = h /^/(g) + o(M- 



We will reduce to this case by means of the Bethe equations in the logarithmic form 



1 ^ ?7rr? ■ 

' fc=l ' 



For a set of variables {^j}i<j<N parametrizing the ground or the quasi-ground state, 
let us introduce the function 



1 N 

9Af(A)=p (A) + — ]Tfl(A-A fc ) 



k=l 



and the set of variables qj = qivi(^j) = 1 < J < A. These new variables have 

a homogeneous distribution. Moreover, qm(A) is an invertible function for any A > 1, 
at least for M > Mo (A) with M large enough. For A > 2, Mo (A) = (and does not 
depend on A). It follows from the estimate 



A' 



q' M {\) = p' (X) - (27T/M) K(\ - X k ) > p' (n/2) - ttK(0) 



k=l 
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(the r.h.s. is positive for A > 2). For 1 < A < 2, the existence of Mo (A) follows from 
the fact that the solution for the ground state does exist and is unique for any M ||[7| , 
and the density p(X) = (l/2ir) lim m ^oo QmW is positive in the thermodynamic limit. 
Note at last that, since gjv/(A + vr) = gM(A) + vr, the function foqjj is 7r-periodic. Hence 
we can represent any sum on j as an integral: 

N N 

iE/( A i) = M^ /fel(5i))i (R2) 

j=i 3=1 



7T 

^ j[dqf(q^( q ))+0(M-^), (B.3) 



2vr 

2 

which we can rewrite as an integral on A: 

N n_ 
T-riZf^) = Y 1 2 d\q' M (\)f(\)+0(M-™). (B.4) 
3=1 ^ J ~% 

One has thus to evaluate the derivative of qu in the thermodynamic limit. We look for 
it in the following form: 

2vr - 

QmW = PoW - J7 E K ( x - = 2 <pW + <PmW), (B.5) 

k=l 

where the correction (fM is to be determined. An integral equation for it can be obtained 
taking /(A) = K(p — A): from ( |B.4| ) one obtains 

■ H 

2 d\K(fi - \)<pmW = -<Pm(h) + 0(M-°°). 



As all the eigenvalues (|4.q ) of the integral operator K are positive, the only function 
satisfying this constraint is up to order 0(M~°°). Hence we proved that 

1 N f~ 

S E/W = ld\ P (A)f(\) + 0(M-°°). (B.6) 

3=1 f 

The proposition 3.2 can be derived directly from this proposition and the analogous 
fact for the homogeneous distribution of momenta. Really for any function g(q) such 
that g'(q) is a periodic function one can easily prove that 

1 N 1 r~ C 

3=1 2 
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(as all further corrections are proportional to the integrals of second and higher deriva- 
tives of g(q) and these integrals are evidently zero). Now one can continue the proof 
as for the first proposition. One should just mention that if g'(X) is a periodic function 
then dg ^ 9 ^ is also a periodic function. 



Appendix C 

In this appendix, we present the calculation of the functions 4>±(X) (3.7) and ip±(\) 
(3.9). 

The functions 4>±(X) can be expressed in the exponential form, 

^(Aa) = exp < ^2 (lnsin(A a - A fc - 5X k ± iQ - lnsin(A a - X k ± iQ) > . 

U=i J (C.l) 

One should note that if all 5X k are replaced by the corresponding 8X k , this expression is 
exactly equal to —1. We will use this property further to calculate the function <^>±(A) 
in the thermodynamic limit. 

Let us at first calculate the "shift functions" 5X k and 5X k . The difference of two 
successive Bethe equations in the logarithmic form ( |2.3D leads to: 

N 

MipoiXj + SXj) - P o(Xj)} + U(Xj + 5Xj - A fe ) - 6{X 3 - X k )\ = 2vr. 

L } k=i } (C.2) 

Representing the differences in the left hand side as Taylor series and using proposition 
l| and Lieb equation ( |3.2[ ), we obtain the following equation for <5A: 



M^—p^-^(\){S\f = l + 0(M-°°). (C.3) 



k=l 



This equation defines uniquely (at order 0(M °°)) the shift function 5X as a functional 

M ' 



of p(X). The solution of this equation can be found as a series on — • 



SX p{X)M i 1 + 2M (p(A)) + 12M 2 (p 2 (A)) + ° (m 3 

We will not use the explicit form of this solution in further calculations. 

The calculation of the function 5X k is more complicated. Taking the difference 
between the Bethe equations in the logarithmic form for the ground state and quasi 
ground state one obtains: 

N 

M{po(Aj + SXj) - poiXj)] + {^j + SXj " Afc - 5~Xk) - d(Xj - Xk)} = 2vr. 

fc=i (C.4) 
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Let us introduce the function z(Xk) = Mp(Xk)SXk- Using the proposition 3J and Lieb 
equation fl3.2|) , one obtains the following equation for this function 

2 

z{\) + J K(X- p)z{p)dp = 1 + (±\ . (C.5) 

7T 

2 

This equation admits a unique solution z(X) = 5 + C^M" 1 ), hence the expression for 
the first order of SXk- 

To obtain the next orders we compare the Bethe equations for the ground and quasi- 
ground state, namely the sums X^fcLi 9{^j — ^k) and YLk=l ^(Aj — Mfc)* The function 9(X) 
is not periodic but its derivative is 7r-periodic, and thus proposition 3.2 can be applied: 

1 / N N \ 

- ^ -/**)- E ^ - Afe ) = £ + 0(M_00) 

\fc=l fc=l / 

where c does not depend on M. As only the first order correction is non-zero, we can 
use merely the first order of 5Xk to calculate c. Finally one obtains 

N N 

£ 0(A, - A*fe) - £ ^ - A fc ) = vr + 0(Ar°°). 
fc=i fc=i 

The Bethe equations for the ground and quasi-ground states can now be written as 

AT 



Mpofjij) + £ 0(jij - n k ) = 2Tmj, (C.6) 

k=l 
N 

MpoiXj) + d ( x j - Ma) = 2 ™j " ^ + 0(M-°°). (C.7) 



fe=l 



Taking the difference of these equations and using once again the Taylor series and the 
Lieb equation we obtain the equation defining 5X as a functional of p(X): 



00 



MV l p ( fc - 1 )(A)(5A) fc = - + 0{M-°°). (C.8) 
z — ' k\ 2 

k=l 

This result is to be compared with flC.3p , It follows that 

SX[p(X)} = SX[2p(X)\ + 0(M~°°). (C.9) 
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To calculate the functions 4>±(X) we use our first remark: replacing 5X k by SX k in 
the expression flC.l|) and developping it as a series on jj leads to 

-1 = exp j^(lnsin(A a - X k - 5X k ± iQ - lnsin(A a - X k ± j , 

\ k — 1 ft- — 1 j 

where j4 n [/?(A)] are homogeneous functionals of p(X) which do not depend on M. The 
first term of the development can be easily computed using proposition 3.1 once again: 

E ( C ° t(A ; ( ^ ±iC) ) = f\ dX cot(X a - X ± iC) + O(M-) = ^ + O(M-). 



fc=i 

Hence it follows that 



A n [p(X)} = 0. 



Thus, due to the relation ( C.9 ), the functions </>±(A) can be represented as the 
following series on -^: 

Since j4 n [p(A)] are homogeneous functionals of p(A), we obtain the very simple result: 

4>±{X) = ±i + 0(M-°°). (CIO) 

The functions ip±(p) can be computed similarly: 

^±(M)=T* + 0(M-°°). (C.ll) 
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